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Abstract 



This is an introduction into the problem of how to set up black hole initial- 
data for the matter-free field equations of General Relativity. The approach is 
semi-pedagogical and addresses a more general audience of astrophysicists 
and students with no specialized training in General Relativity beyond that 
of an introductory lecture. ' 



1 Introduction and motivation 

In my lecture I will try to explain how scattering and merging processes between 
black holes can be described analytically in general relativity(GR). This is a vast 
subject and I will focus attention to the basic issues, rather than trying to explain 
analytical details of approximation schemes etc. I will also not discuss numerical 
aspects, which are beyond my competence, and which would anyway require a 
separate lecture. I will address the following main topics: 

1) a first step beyond Newtonian gravity, 

2) constrained evolutionary structure of Einstein's equations, 

3) the 3-1-1- split and the Cauchy initial- value problem, 

4) black hole data, 

5) problems and recent developments. 



'This is the written version of a lecture delivered at the school "The Galactic Black Hole", held 
between August 26.-3 1 . in 2001 at the Physics Center of the German Physical Society in Bad Honnef 
(Germany). It is published in: Heino Falcke and Friedrich W Hehl (editors) The Galactic Black Hole 
(Lectures on General Relativity and Astrophysics), lOP Publishing (Bristol) 2003. 
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with emphasis on the fourth entry. However, I will also spend some time to explain 
some of the specialties of GR, like the absence of a point-particle concept and 
the non-trivial linkage between the field equations and the equations of motion for 
matter. These points should definitely be appreciated before one goes on to discuss 
black holes, which are solutions to the vacuum Einstein equations representing 
extended objects. 

The following points seem to me the main motivations for studying the problem 
of black hole collision: 

• Coalescing black holes are regarded as promising sources for the detection 
of gravitational waves by earth-based instruments. 

• Close encounters of black holes provide physically relevant situations for the 
investigation of the strong-field regime of general relativity. 

• The dynamics of simple black hole configurations is regarded as ideal testbed 
for numerical relativity. 

My conventions are as follows: space-time is a manifold M with Lorentzian metric 
g of signature (—,+,+,+). Greek indices range from to 3, latin indices from 1 to 
3 unless stated otherwise. The covariant derivative is denoted V^^, ordinary partial 
derivatives by dfj, or sometimes simply by a lower-case , /x. The relation := (=:) 
defines the left (right) hand side. The gravitational constant in GR is k = SvrG/c^, 
where G is Newton's constant and c the velocity of light. A symbol hke (9(e") 
stands collectively for terms falling off at least as fast as e". 

2 A first step beyond Newtonian gravity 

It can hardly be overstressed how useful the concept of a mass-point is in New- 
tonian mechanics and gravity. It allows to pointwise probe the gravitational field 
and to reduce the dynamical problem to the mathematical problem of finding so- 
lutions to a system of finitely many ordinary differential equations. To be sure, 
just postulating the existence of mass points is not sufficient. To be consistent with 
known laws of physics one must eventually understand the point mass as an ideal- 
ization of a highly locahzed mass distribution which obeys known field-theoretic 
laws, such that in the situations at hand most of the field degrees of freedom effec- 
tively decouple from the dynamical laws for those collective degrees of freedom 
one is interested in, hke e.g. the centre-of-mass. In Newtonian gravity this usually 
requires clever approximation schemes but is not considered to be a problem of 
fundamental nature. Although this is true for the specific hnear theory of Newto- 
nian gravity, this needs not be so for comparably simple generalizations, as will 
become clear below. 
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In GR the situation is markedly different. A concentration of more than one 
Schwarzschild mass in a region of radius less than the Schwarzschild radius will 
lead to a black hole whose behaviour away from the stationary state usually can- 
not be well described by finitely many degrees of freedom. It shakes and vibrates, 
thereby radiating off energy and angular momentum in form of gravitational radi- 
ation. Moreover, it is an extended object and cannot be unambiguously ascribed 
an (absolute or relative) position or individual mass. Hence the problem of mo- 
tion, and therefore the problem of scattering of black holes, cannot be expected to 
merely consist of corrections to Newtonian scattering problems. Rather, the whole 
kinematic and dynamical setup will be different where many of the established 
concepts of Newtonian physics need to be replaced or at least adapted, very often 
in a somewhat ambiguous way. Among those are mass, distance, and kinetic en- 
ergy. For example, one may try to solve the following straightforward sounding 
problem in GR, whose solution one might think has been given long ago. Given 
are two unspinning black holes, momentarily at rest, with equal individual mass m, 
mutual distance i, and no initial gravitational radiation around. What is the amount 
of energy released via gravitational radiation during the dynamical infall? In such 
situation we can usually make sense of the notions of 'spin' (hence unspinning) and 
'mass'; but ambiguities generally exist in defining 'distance' and, most important 
of all, 'initial gravitational radiation'. Such difficulties persist over and above the 
ubiquitous analytical and/or numerical problems which are currently under attack 
by many research groups. 

To those who are not so familiar with GR and like to see Newtonian analogies, 
I wish to mention that there is a way to consistently model some of the non-linear 
features of Einstein's equations in a Newtonian context, which shares the property 
that it does not allow for point masses. I will briefly describe this model since it 
does not seem to be widely known. 

First recall the field equation in Newtonian gravity, which allows to determine 
the gravitational potential (whose negative gradient, -Vcp, is the gravitational 
field) from the mass density p (the 'source' of the gravitational field): 

A^ = iTTGp. (1) 

Now suppose one imposes the following principle for a modification of Q: all 
energies, including the self-energy of the gravitational field, act as source for the 
gravitational field. In order to convert an energy density s into a mass density p, 
we adopt the relation e = p(? from special relativity (the equation we will arrive 
at can easily be made Lorentz invariant by adding appropriate time-derivatives). 
The question then is whether one can modify the source term of Q such that p 
P+PgTs.^ with pgrav := Egrav/c^, whcrc Sg^av is the energy density of the gravitational 
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field as predicted by this very same equation (condition of self-consistency). It 
turns out that there is indeed a unique such modification, which reads 



47rG 



P + 




(V<A/0)2 



(2) 



It is shown in the Appendix that this equation indeed satisfies the 'energy principle' 
as just stated. (For more information and a proof of uniqueness, see Oil .) The 
gravitational potential is now required to be always positive, tending to the value 
at spatial infinity (rather than zero as for (Q). The second term on the right hand 
side of ^ corresponds to the energy density of the gravitational field. Unlike the 
energy density following from © (which is — g^l^'^l^-' '^"^ positive definite. 
This does not contradict attractivity of gravity for the following reason: the rest- 
energy density of a piece of matter is in this theory not given by pc^, but by pcj), that 
is, it depends on the value of the gravitational potential at the location of matter. 
The same piece of matter located at lower gravitational potential has less energy 
than at higher potential values. In GR this is called the universal redshift effect. 
Here, as in GR, the active gravitational-mass also suffers from this redshift, as is 
immediate from the first term on the right hand side of where p does not enter 
alone, as in Q, but is multiplied with the gravitational potential 0. With respect to 
these features our modification ^ of Newtonian gravity mimics GR quite well. 

We mention in passing that ^ can be 'linearized' by introducing the dimen- 
sionless field ijj := y/c^/c, in terms of which ^ reads 



combinations of solutions are again solutions whose coefficients add up to one. For 
p > it also follows that solutions to ^ can never assume negative values, since 
otherwise the function ijj must have a negative minimum (because of the positive 
boundary values) and therefore non-negative second derivatives there. But then ^ 
cannot be satisfied at the minimum, hence if) must be nonnegative everywhere. This 
implies that solutions of ^ are also non-negative. To be sure, for mathematical 
purposes ^ is easier to use than but note that cj) and not ijj is the physical 
gravitational potential. 

We now show how these non-linear features render impossible the notion of 
point mass, and even induce a certain black-hole behaviour on their solutions. Let 
us be interested in static, spherically symmetric solutions to ^ with source p, 
which is zero for r > i? and constant for r < R. We need to distinguish two 
notions of mass. One mass just counts the amount of 'stuff located within r < R. 




(3) 



The boundary conditions are now ip{r 



1. Hence only those linear 
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You may call it 'bare mass' or 'baryonic mass', since for ordinary matter it is 
proportional to the baryon number. We denote it by Mb- It is simply given by 



Mb := / d?xp. (4) 

J space 

The other mass is the 'gravitational mass', which is measured by the amount of 
flux of the gravitational field to 'infinity', that is, through the surface of a sphere 
whose radius tends to infinity. We call this mass Mq. It is given by 



Mg :=-^lim/ {-V<t)-n)da, (5) 

4vrG r->oo /_5.2(^) 



->oo 



where r = \x\, x/r = n, S'^{r) is the 2-sphere or radius r and da is its surface 
element. Mq should be identified with the total inertial mass of the system, in full 
analogy to the ADM-mass in GR (see eq. d35l below). Hence Mqc^ is the total 
energy of the system, with gravitational binding energy also taken into account. 
The masses Mb and Mq can dimensionally be turned into radii by writing Rb '■= 
GMbI(? and Rq := GMg/c^, and further turned into dimensionless quantities 
via rescaling with R, the radius of our homogeneous star. We write x := Rb/ R 
and y := Rg/R- 

For each pair of values for the two parameters AIb and R there is a unique 
homogeneous-star-solution to whose simple analytical form needs not inter- 
est us here (see |21|). Using it we can calculate Mq, whose dependence on the 
parameters is best expressed in terms of the dimensionless quantities x and y: 



tanh(73^) ^ 

y = f^^^ = '[' — Ti^j- 

The function / maps the interval [0, oo] monotonically to [0, 2]. This implies the 
following inequality 

Mg < 2Rc^/G (7) 

which says that the gravitational mass of the star is bounded by a purely geometric 
quantity. It corresponds to the statement in GR that the star's radius must be big- 
ger than its Schwarzschild radius, which in isotropic coordinates is indeed given 
by Rs = GMgIIc?. It can be proven ED that the bound Q still exists for 
non-homogeneous spherically symmetric stars, so that the somewhat unphysical 
homogeneity assumption can be lifted. The physical reason for this inequality is 
the 'redshift', i.e. the fact, that the same bare mass at lower gravitational poten- 
tial produces less gravitational mass. Hence adding more and more bare mass into 
the same volume pushes the potential closer and closer to zero (recall that <^ is 
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always positive) so that the added mass becomes less and less effective in gener- 
ating gravitational fields. The inequality then expresses the mathematical fact that 
this 'redshifting' is sufficiently effective so as to give finite upper bounds to the 
gravitational mass, even for unbounded amounts of bare mass. 

The energy balance can also be nicely exhibited. Integrating the matter en- 
ergy density (pp and the energy density of the gravitational field, -^^{V (j) / cj))^ , we 
obtain 

^.atte. = Mbc^ " ^ + OiRl/R")^ , (8) 

i^flcM = Mbc" (^^ + 0{Rl/R')^ , (9) 

i^tota, = Mbc" (^l-^ + 0{Rl/R^)^ =Mgc\ (10) 

Note that the term —SMbc^Rb/^R in •QOJ is just the Newtonian binding energy. 
At this point it is instructive to verify the remarks we made above about the posi- 
tivity of the gravitational energy. Shrinking a mass distribution enhances the field 
energy, but diminishes the matter energy twice as fast, so that the overall energy 
is also diminished, as it must due to the attractivity of gravity. But here this is 
achieved with all energies involved being positive, unlike in Newtonian gravity. 
Note that the total energy, Mq, cannot become negative (since (p cannot become 
negative, as already shown). Hence one also cannot extract an infinite amount of 
energy by unlimited compression, as it is possible in Newtonian gravity. This is 
the analogue in our model theory to the positive mass theorem in GR. 

We conclude by making the point announced above, namely that the inequality 
(Q shows that point objects of finite gravitational mass do not exist in the the- 
ory based upon (|2li; mass implies extension! Taken together with the lesson from 
special relativity, that extended rigid bodies also do not exist (since the speed of 
elastic waves is less than c), we arrive at the conclusion that the dynamical prob- 
lem of gravitating bodies and their interaction is fundamentally of field-theoretic 
(rather than point-mechanical) nature. Its proper realization is GR to which we 
now turn. 



3 Constrained evolutionary structure of Einstein's equa- 
tions 

In GR the basic field is the spacetime metric g^i^, which comprises the gravitational 
and inertial properties of spacetime. It defines what inertial motion is, namely a 
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geodesic 

+ r^^x^i'' = (11) 
with respect to the Levi-Civita connection 

r^i. := h^''{-9t^u,a + gai,,u + gua,)^) ■ (12) 

(Since inertial motion is 'force-free' by definition, you may rightly ask whether 
it is correct to call gravity a 'force'.) The gravitational field g^,^ is linked to the 
matter content of spacetime, represented in form of the energy momentum tensor 
T^jy, by Einstein's equations 

G^i/ '.= R^u — 29fj.i'R — i^T^y. (13) 

Due to the gauge-invariance with respect to general differentiable point transforma- 
tions (i.e. diffeomorphisms) of spacetime one has the identities (as a consequence 
of Noether's 2nd theorem) 

V^G^"' = . (14) 

Being 'identities' they hold for any G^'^ , independent of any field equation. With 
respect to some coordinate system = {x^ , • • • x^) we can expand (I14t in terms 
of ordinary derivatives. Preferring the coordinate x", this reads 

d^G^'^ = -dkG^'' - T^^y^G^'' - T^'^xG^^- (15) 

Since G^''^ contains no higher derivatives of g^j^y than the second, the right-hand 
side of this equation also contains only 2nd x'^-derivatives. Hence (fTSt implies 
that the four components G^^ only involve first x'^-derivatives. Now choose x° as 
time coordinate. The four (0, zv) -components of (fT3t then do not involve second 
time derivatives, unlike the space-space components Hence the time-time 

and time-space components are constraints, that is, equations that constrain the 
allowed choices of initial data, rather than evolving them. 

That is not an unfamiliar situation, which similarly occurs for Maxwell's equa- 
tions in electromagnetism (EM). Let us recall this analogy. We consider the four- 
dimensional form of Maxwell's equations in terms of the vector potential A^, 
whose antisymmetric derivative is the field-tensor F^y := d^Ay — dyA^, compris- 
ing the electric {Ei = F^i) and magnetic {Bi = —Fj^, ijk cyclic) field. Maxwell's 
equations are 

E-^ -d^F'^-' -^f = 0. (16) 
Due to its antisymmetry, the field tensor obviously obeys the identity 

d^dyFf''' = 0, (17) 
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which here is the analogue of (fT4b . an identity involving third derivatives in the 
field variables. Using ( 1171 in the divergence of ( fT^ yields 



d,E^ = -^d,j\ (18) 

which shows that Maxwell's equations imply charge conservation as integrabil- 
ity condition. Let us interpret the role of charge conservation in the initial value 
problem. Decomposing (flTl into space and time derivatives gives 

dod.F^" = -dkd.F'^'' . (19) 

Again the right hand side involves only second time derivatives implying that the 0- 
component of (IT^ involves no second time derivatives. Hence the time component 
of (fT6t is merely a constraint on the initial data; clearly it is just Gauss' law V • 
E — 47rp = 0. Its change under time evolution according to Maxwell's equations 
is 

ao^° = d.E'^-dkE^ 

= -^d,f-dkE\ (20) 

where we used the identity ( fTSl in the second step. Suppose now that on the initial 
surface of constant we put an electromagnetic field which satisfies the con- 
straint, E^ = 0, and which we evolve according to E''^ = (implying d^E^ = Q 
on that initial surface). Then d20l shows that charge conservation is the necessary 
and sufficient condition for the evolution to preserve the constraint. 

Let us return to GR now, where the overall situation is entirely analogous. Now 
we have four constraints 

^Ou .-qOu _ ^j,Ou ^ Q (-21) 

and six evolution equations, which we write 

E'^ := G'^ - kT^ = Q . (22) 
The identity (fT?t now implies 

V;,^^^ = -kV^T^'' , (23) 



which parallels (TTSl . Here, too, the time derivative of the constraints is easily 
calculated: 

dnE^'' = V.Ef"' -di.E^'' -T^a^E^'' -n^E^^ 
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using (l23l in the last step. Now consider again the evolution of initial data from 
a surface of constant x^. If they initially satisfy the constraints and are evolved 
via E^^ = (hence all spatial derivatives of E'^'^ vanish initially) they continue 
to satisfy the constraints if and only if the energy momentum tensor of the matter 
satisfies 



Hence we see that the 'covariant conservation' of energy-momentum, expressed by 
(l25l . plays the same role in GR as charge conservation plays in EM. This means 
that you cannot just prescribe the motion of matter and then use Einstein's equa- 
tions to calculate the gravitational field produced by that source. You have to move 
the matter in such a way that it satisfies (l25l . But note that at this point there is a 
crucial mathematical difference to charge conservation in EM: charge conservation 
is a condition on the source only, it does not involve the electromagnetic field. This 
means that you know a priori what to do in order not to violate charge conserva- 
tion. On the other hand, (l25l involves the source and the gravitational field. The 
latter enters through the covariant derivatives which involve the metric g^i^ through 
the connection coefficients (IT2l . Hence here (l25l cannot be solved a priori by suit- 
ably restricting the motion of the source. Rather we have a consistency condition 
which mutually links the problem of motion for the sources and the problem of 
field determination. It is this difference which makes the problem of motion in 
GR exceedingly difficult. (A brief and lucid presentation of this problem, drawing 
attention to its relevance in calculating the generation of gravitational radiation by 
self-gravitating systems, was given in ifTSll . A broader summary, including modem 
developments is Hl^ ). 

For example, for pressureless dust represented by T^^'^ = p(? U^V^ , where 
p is the local rest-mass density and C/^ is the vector field of four-velocities of 
the continuously dispersed individual dust grains, (l25l is equivalent to the two 
equations 



The first states the conservation of rest-mass. The second is equivalent to the state- 
ment that the vector field C/^ is geodesic, which means that its integral lines (the 
worldlines of the dust grains) are geodesic curves (IIU with respect to the met- 
ric g^y. Hence we see that in this case the motion of matter is fully determined 
by d25l . i.e. by Einstein's equations, which imply (l25l as integrability condition. 
This clearly demonstrates how the problem of motion is inseparably linked with 
the problem of field determination, and that these problems can only be solved 
simultaneously. The methods used today use clever approximation schemes. For 




(25) 



V^ipU^') = 0, 
U^VyW = 0. 



(26) 
(27) 
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example, one can make use of the fact that there is a difference of one power in k 
between the field equations and their integrability condition. Hence, in an approxi- 
mation in K it is consistent for the nth order approximation of the field equations to 
have the integrability conditions (equations of motions) satisfied to n — 1st order. 

Clearly the problem just discussed does not arise for the matter-free Einstein 
equations for which T^,^ = 0. Now recall that black holes too are described by 
the matter free equations. Hence the mathematical problem just described does not 
occur in the discussion of their dynamics. In this aspect the discussion of black 
hole scattering is considerably easier than e.g. that of neutron stars. 

4 The 3+1- split and the Cauchy initial- value problem 

We saw that the ten Einstein equations decompose into two sets of four and six 
equations respectively, four constraints which the initial data have to satisfy, and 
six equations driving the evolution. As a consequence there will be four dynam- 
ically undetermined components among the ten components of the gravitational 
field gfj_^. The task is to parametrise the g^,^ in such a way that four dynamically 
undetermined functions can be cleanly separated from the other six. One way to 
achieve this is via the splitting of spacetime into space and time (see |22| for a 
more detailed discussion). The four dynamically undetermined quantities will be 
the famous lapse (one function a) and shift (three functions /?*). The dynamically 
determined quantity is the Riemannian metric hij on the spatial 3-manifolds of 
constant time. These together parametrise g^,^ as follows: 

= _c^2 (^^0)2 ^ ^^^t ^ p^dx^)^dx'' + [3^ dx^) . (28) 

The physical interpretation of a and /3* is as follows: think of spacetime as the 
history of space. Each 'moment' of time, x^ = t, corresponds to an entire 3- 
dimensional slice E^. Obviously there is plenty of freedom how to 'waft' space 
through spacetime. This freedom corresponds precisely to the freedom to choose 
the 1-1-3 functions a and For one thing, you may freely specify how far for each 
parameter step dt you push space in a perpendicular direction forward in time. This 
is controlled by a, which is just the ratio ds / dt of the proper perpendicular distance 
between the hypersurfaces and Tit+dt- This speed may be chosen in a space and 
time dependent fashion, which makes a a function on spacetime. Second, given 
a point with coordinates on St. Going from x* in a perpendicular direction 
you meet T^t+dt in a point with coordinates + dxi, where dx* can be chosen 
at will. This freedom of moving the coordinate system around while evolving is 
captured by one writes dx* = f3^dt. Clearly this moving around of the spatial 
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coordinates can also be made in a space and time dependent fashion, so that the /3* 
are functions of spacetime, too. 

Let n^^ be the vector field in spacetime which is normal to the spatial sections of 
constant time. It is given by n = -{d/dx^ — jS^d/dx^), as one may readily verify 
using ( EHl (you have to check that n is normalized and satisfies g{n, d/dx^) = 0). 
We define the extrinsic curvature, Kij, to be one-half the Lie derivative of the 
spatial metric in the direction of the normal: 



where D is the spatial covariant derivative with respect to the metric hij . As usual, 
a round bracket around indices denotes their symmetrization. Note that, by defini- 
tion, Kij is symmetric. Finally we denote the Ricci scalar of h^j by R^'^^K 

We can now write down the four constraints of the vacuum Einstein equations 
in terms of these variables: 



(1301) and (I31t are referred to as Hamiltonian constraint and momentum constraint 
respectively. The six evolution equations of second order in the time derivative 
can now be written as twelve equations of first order. Six of them are just ( l29t . 
read as equation that relates the time derivative dhij/dx^ to the 'canonical data' 
{hij, Kij). The other six equations, whose explicit form needs not concern us here 
(see e.g. |22|), express the time derivative of Kij in terms of the canonical data. 
Both sets of evolution equations contain on their right hand sides the lapse and 
shift functions, whose evolution is not determined but must be specified by hand. 
This specification is a choice of gauge, without which one cannot determine the 
evolution of the physical variables {hij, Kij). 

The initial-data problem takes now the following form 

I. Choose a topological 3-manifold S. 

II. Find on S a Riemannian metric hij and a symmetric tensor field Kij which 
satisfy the constraints d30t and dSTT i. 

III. Choose a lapse function a and a shift vector field /5\ both as functions of 
space and time, possibly according to some convenient prescription (e.g. 
singularity avoiding gauges, like maximal slicing). 




(29) 



= G{n,n) = ^{R^'^ + K'^Kij-{Kl)^) 
= G{n, d/dx^) = Di{m - d'jK^) . 



(30) 
(31) 
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IV. Evolve initial data with these choices of a and /?' according to the twelve 
equations of first order. By consistency of Einstein's equations the con- 
straints will be preserved during this evolution, independent of the choices 
for a and 

The backbone of this setup is a mathematical theorem, which states that for any set 
of initial data, taken from a suitable function space, there is, up to diffeomoiphism, 
a unique maximal Einstein spacetime developing from these data L2J- 

5 Black hole data 
5.1 Horizons 

By black hole data we understand vacuum data which contain apparent horizons. 
The informal definition of an apparent horizon is that it is the boundary of a trapped 
region, which means that its orthogonal outgoing null rays must have zero diver- 
gence. (Inside the trapped region they converge for any 2-surface, by definition of 
'trapped region'.) The Penrose-Hawking singularity theorems state that the exis- 
tence of an apparent horizon implies that the evolving spacetime will be singular 
(assuming the strong energy-condition). Given also the condition that singularities 
cannot be seen by observers far off, a condition usually called cosmic censorship, 
one infers the existence of an event horizon and hence a black hole. One can then 
show that the intersection of the event horizon with the spatial hypersurface lies 
on or outside the apparent horizon (for stationary spacetimes they coincide). The 
reason why one does not deal with event horizons directly is that you cannot tell 
whether there exists one by just looking at initial data. In principle you would have 
to evolve them to the infinite future, which is beyond our abilities in general. In 
contrast, apparent horizons can be recognized once the data on an initial slice are 
given. The formal definition of an apparent horizon is the following: given initial 
data (S, hij,Kij) and an embedded 2-surface a C S with outward pointing nor- 
mal I/*. (T is an apparent horizon if and only if the following relation between Kij, 
the extrinsic curvature of S in spacetime, and kij, the extrinsic curvature of a in S, 
is satisfied: 

qiJkij = -q'^Kij , (32) 

where qij := hij — ViVj is just the induced Riemannian metric on a, so that (l32t 
simply says that the restriction of Kij to the tangent space of a has opposite trace 
to kij. (The minus sign on the rhs of d32t signifies future apparent horizon cor- 
responding to a black hole which has a future event horizon. A plus sign would 
signify a past apparent horizon corresponding to a 'white hole' with past event 
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horizon.) This means that once we have the data (S, hij,Kij) we can in princi- 
ple find all 2-surfaces o" C S for which &7i holds and therefore find all apparent 
horizons. 

5.2 Poincare charges 

By Poincare charges we shall understand quantities like mass, linear momentum, 
and angular momentum. In GR they are associated to an asymptotic Poincare 
symmetry (see |4|), provided that the data {T,,hij,Kij) are asymptotically flat in 
a suitable sense, which we now explain. Topologically asymptotic flatness means 
that the non-trivial 'topological features' of S should all reside in a bounded region 
and not 'pile up' at infinity. More formally this is expressed by saying that there is 
a bounded region B C Tj such that S — i? (the complement of B) consists of a finite 
number of disjoint pieces, each of which looks topologically like the complement 
of a ball in M^. These asymptotic pieces are also called the ends of the manifold 
S. Next comes the geometric restriction imposed by the condition of asymptotic 
flatness. It says that for each end there is an asymptotically euclidean coordinate- 
system {x^, , x^} in which the fields {hij, Kij) have the following fall-off for 
r ^ oo (r = [x^Y + (3^^)^ + {x'^Y, = x^ /r) 

h,,{x^) = 5,, + ^^ + 0(r-i-^), (33) 
K^.ix'^) = Mj!) +0(^-2-.). (34) 

Moreover, in order to have convergent expressions for physically relevant quanti- 
ties, like e.g. angular momentum (see below), the field Sij must be an even func- 
tion of its argument, i.e. Sij{—n^) = Sij{n^), and tij must be an odd function, i.e. 

Uj{-n^) = -Uj{n^). 

Under these conditions each end can be assigned mass, momentum, and angu- 
lar momentum, which are conserved during time evolution. They may be computed 
by integrals over 2-spheres in the limit the spheres are pushed to larger and larger 
radii into the asymptotically flat region of that end. These so-called ADM-Integrals 
(first considered by Arnowitt, Deser, and Misner in 1 1 1) are given by the following 
expressions, which we give in 'geometric' units (meaning that in order to get them 
in standard units one has to multiply the mass expression given below by 1/k and 
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the linear and angular momentum by c/k): 



M 



lim / S^Hdihjk — dkhij)n^ da 



(35) 




Js^{r) 



(36) 



(37) 



5.3 Maximal and time-symmetric data 

The constraints ( OUIOTI are too complicated to be solved in general. Further condi- 
tions are usually imposed to reduce the complexity of the problem: data {hij , Kij ) 
are called maximal if = h^^ Kij = 0. The name derives from the fact that 
Kl = is the necessary and sufficient condition for a hypersurface to have station- 
ary volume to first order with respect to deformations in the ambient spacetime. 
Even though stationarity does generally not imply extremality one calls such hy- 
persurfaces maximal. Note also that since spacetime is a Lorentzian manifold, 
extremal spacelike hypersurfaces will be of maximal rather than minimal volume. 
In contrast, in Riemannian manifolds one would speak of minimal surfaces. 

A much stronger condition is to impose Kij = 0, which as seen from Mb\ and 
(l37l implies that all momenta and angular momenta vanish. Only the mass is now 
allowed to be non-zero. Such data are called time-symmetric since for them hij 
is momentarily static as seen from ( 129 1 . This implies that the evolution of such 
data into the future and into the past will coincide so that the developed space- 
time will have a time-reversal symmetry which pointwise fixes the initial surface 
where Kij = 0. This surface is therefore also called the moment of time-symmetry . 
Time-symmetric data can still represent configurations of any number of black 
holes without angular momenta which are momentarily at rest. Note also that for 
time-symmetric data the condition (l32t for an apparent horizon is equivalent to 
the tracelessness of the extrinsic curvature of a. Hence for time symmetric data 
apparent horizons are minimal surfaces. 

We add one more general comment concerning submanifolds. A vanishing ex- 
trinsic curvature is equivalent to the property that each geodesic of the ambient 
space, which starts on, and tangent to, the submanifold, will always run entirely 
inside the submanifold. Therefore, submanifolds with vanishing extrinsic curva- 
ture are called totally geodesic. Now, if the ambient space allows for an isometry 
(symmetry of the metric), whose fixed-point set is the submanifold in question, 
like for the time-reversal transformation just discussed, the submanifold must nec- 
essarily be totally geodesic. To see this, consider a geodesic of the ambient space 
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which starts on, and tangent to, the submanifold. Assume that this geodesic even- 
tually leaves the submanifold. Then its image under the isometry would again be a 
geodesic (since isometrics always map geodesies to geodesies) which is different 
from the one we started from. But this is impossible since they share the same 
initial conditions which are known to determine the geodesic uniquely. Hence 
the geodesic cannot leave the submanifold, which proves the claim. We will later 
have more opportunities to identify totally geodesic submanifolds — namely appar- 
ent horizons — by their property of being fixed-point sets of isometrics. 

5.4 Solution strategy for maximal data 

Possibly the most popular approach to solving the constraints is the conformal 
technique due to York, Lichnerowicz and others (see l37l for a review). The basic 
idea is to regard the Hamiltonian constraint d30l as equation for the conformal 
factor of the metric hij w& freely specify the complementary information, called 
the conformal equivalence-class of hij. More concretely, this works as follows: 

I. Choose unphysical ('hatted') quantities (hij,Kij), where hij is a Rieman- 
nian metric on E and Kij is symmetric, trace and divergence free: 

h'^Kij = 0, D'kij = 0, (38) 

where D is the covariant derivative with respect to hij. 

II. Solve the (quasilinear elliptic) equation for a positive, real valued function 
$ with boundary condition <I>(r — > oo) 1, where A = h^^ DiDj: 

A<^ + Ik'^ kij'^-'^ = (39) 

III. Using the solution of i39l . define physical ('unhatted') quantities by 

hij = ^%j, (40) 
Kij = ^-^ki,. (41) 

These will satisfy the constraints ( 1301 IJil l.^ 

5.5 Explicit time symmetric data 

Before we say a little more about maximal data we wish to present some of the most 
popular examples for time symmetric data some of which are also extensively used 
in numerical simulations. Hopefully these examples let you gain some intuition 
in the geometries and topologies involved and also let you anticipate the richness 
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that a variable space-structure gives to the solution space of one of the simplest 
equations in physics: the Laplace equation. 

Restricting the solution strategy, outlined above, to the time-symmetric case 
one first observes that for Kij = one has Kij = 0. The momentum constraint 
(Hn is automatically satisfied and all that remains is equation ( B^ . which now 
simply becomes the Laplace equation for the single scalar function ^> on the Rie- 
mannian manifold (S, /ijj). 

We now make a further simplifying assumption, namely that hij is in fact the 
flat metric. This will restrict our solution hij to a confomially flat geometry. It is 
not obvious how severe the loss of physically interesting solutions is by restricting 
to conformally flat metrics. But we will see that the latter already contain many 
interesting and relevant examples. 

So let us solve Laplace's equation in flat space! Remember that $ must be 
positive and approach 1 at spatial infinity (asymptotic flatness). We cannot take 
E = since the only solution to the Laplace equation in which asymptotically 
approaches 1 is identically 1. We must allow $ to blow up at some points, which 
we can then remove from the manifold. In this way we let the solution tell us 
what topology to choose in order to have an everywhere regular solution. You 
might think that just removing singular points would be rather cheating, since the 
resulting manifold may turn out to be incomplete, that is, can be hit by a curve after 
finite proper length (you can go 'there'), even though <I> and hence the physical 
metric blows up at this point. If this were the case one definitely had to say what a 
solution on the completion would be. But, as a matter of fact, this cannot happen 
and the punctured space will turn out to be complete in the physical metric. 



5.5.1 One black hole 

The simplest solution with one puncture (at r = 0) is just 

^{r,e,if) = l + -, (42) 
r 

where a is a constant which we soon interpret and which must be positive in order 
for $ to be positive everywhere. We cannot have other multipole contributions 
since they inevitably would force <1> to be negative somewhere. What is the geom- 
etry of this solution? The physical metric is 

ds"^ = (^ + ^'^ {dr^ + dO'^ + W 9 d(p^) (43) 

which is easily checked to be invariant under the inversion-transformation on the 
sphere r = a: 

r^—, e^e, if^ip. (44) 
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Figure 1 : One black hole 



This means that the region r > a just looks Uke the region r < a and that the 
sphere r = a has the smallest area among all spheres of constant radius. It is a 
minimal surface, in fact even a totally geodesic submanifold, since it is the fixed 
point set of the isometiy d44b . Hence it is an apparent horizon, whose area follows 
from (l43l : 

A = 167r(2a)2 . (45) 

Our manifold thus corresponds to a black hole. Its mass can easily be computed 
from (I35t : one finds m = 2a. This manifold has two ends, one for r ^ oo and 
one for r ^ 0. They have the same geometry and hence the same ADM mass, as it 
must be the case since individual and total mass clearly coincide for a single hole. 

The data just written down correspond to the 'middle' slice right across the 
Kruskal (maximally extended Schwarzschild) manifold. Also, (l43b is just the spa- 
tial part of the Schwarzschild metric in isotropic coordinates. Hence we know its 
entire future development in analytic form. Already for two holes this is not the 
case anymore. Even the simplest two body problem — head-on collision — has not 
been solved analytically in GR. 

5.5.2 Two black holes 

There is an obvious generalization of (l42l by allowing two 'monopoles' of strength 
ai and 02 at the punctures x = xi and x = X2 respectively. The 3-metric then reads 

ds^=(l + + V idr' + dd^ + sin2 9 d^'). (46) 

\ — Xl| \X — X2\ J 

The manifold has now three asymptotically flat ends, one for \x\ -^00, where the 
overall ADM mass M is measured, and one each for \x — xi^2\ 0. To see the 
latter, it is best to write the metric (l46l in spherical polar coordinates (ri ,9i,ipi) 
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Figure 2: Two black holes well separated 



centered at xi, and then introduce the inverted radial coordinate given by ri = 
af/ri- In the limit fi ^ oo the metric then takes the form 

(47) 

where ri2 = |xi — X2|- This looks just Uke a one-hole metric (l43l . Hence, if 
the black holes are well separated (compared to their size), the two-hole geometry 
looks like that depicted in figure |2l By comparison with the one-hole metric we 
can immediately write down the ADM masses corresponding to the three ends 
n,2 — ^ oo respectively: 

M = 2(01+02), mi,2 = 201,2(1 + Xi,2), where ^1,2 = — • (48) 

Momenta and angular momenta clearly vanish (moment of time symmetry). StiU 
assuming well separated holes, i.e. Xi = cii/r 12 <C 1, we can calculate the binding 
energy A£' = M — mi — m2 as function of the masses nii and ri2 and get 

= -™ ( 1 - + Oiimi,2/n2f)) . (49) 

ri2 V 2ri2 y 

The leading order is just the Newtonian expression for the binding energy of two 
point-particles with masses mi 2 at distance ri2. But there are corrections to this 
Newtonian form which tend to diminish the Newtonian value. Note also that (l49l 
is still not in a good form since ri2 is not an invariantly defined geometric distance 
measure. As such one might use the length i of the shortest geodesic joining the 
two apparent horizons 5*1 and ^2. Unfortunately these horizons are not easy to 
locate analytically and hence no closed form of i{rni,m2,ri2) exists which could 
be inverted to eliminate ri2 in favour of £. 

Due to the difficulty to analytically locate the two apparent horizons we also 
cannot write down an analytic expression for their area. But we can give upper and 
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Figure 3: Two black holes after merging 



lower bounds as follows: 

16TT{2aif <Ai< 167r[(2ai(l + Xi)? = l^irmf . (50) 

The lower bound simply follows from the fact that the two-hole metric d?^ . if 
written down in terms of spherical polar coordinates about any of its punctures, 
equals the one hole metric (H^ plus a positive definite correction. The upper bound 
follows from the so-called Penrose inequality in Riemannian Geometry (proven in 
||27l), which directly states that IGvrm^ > A for each asymptotically flat end, 
where m is the mass according to d35l and A is the area of the outermost (as seen 
from that end) minimal surface. 

If the two holes approach each other to a distance comparable to the sizes of 
the holes the geometry changes in an essential way. This is shown in figure |3] The 
most important new feature is that new minimal surfaces form, in fact two 1 10|, 
which both enclose the two holes. The outermost of these, as seen from the upper 
end, denoted by S3 in figure |3l corresponds to the apparent horizon of the newly 
formed 'compound' black hole which contains the two old ones. For two black 
holes of equal mass, i.e. ai = 02 = o, this happens approximately for a parameter 
ratio of a/ri2 = 0.65, which in an approximate numerical translation into the ratio 
of individual hole mass to geodesic separation reads m/£ 0.26. 

5.5.3 More than two black holes 

The method can be generalized in a straightforward manner to any number n of 
black holes with parameters (oj, Xj), i = 1, • • • n, for the punctures. The manifold 
S has now n + 1 ends, one for |x| — > 00 and one for each x ^ Xi. The expressions 
for the metric and masses are then given by the obvious generahzations of (l46t and 
(l48l respectively. 
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5.5.4 Energy bounds from Hawking 's area law 



Loosely speaking, Hawking's area law states that the surface of a black hole cannot 
decrease with time. (See 1^ for a simple and complete outline of the traditional 
and technically slightly restricted version and 19] for the technically most complete 
proof known today.) Let us briefly explain this statement. If S is a Cauchy surface 
(a spacelike hypersurface in spacetime) and Ti the event horizon (a lightlike hy- 
persurface in spacetime), the two intersect in a number of components (spacelike 
2-manifolds), each of which we assume to be a two-sphere. Each such two-sphere 
is called the surface of a black hole at time S. Let us pick one of them and call it 
B. Consider next a second Cauchy surface S' which lies to the future of S. The 
outgoing null rays of B intersect S' in a surface B' , and the statement is now that 
the area of B' is larger than or equal to the area of B (to prove this one must assume 
the strong energy condition). Note that we deliberately left open the possibility that 
B' might be a proper subset of a black hole surface at time S', namely in case the 
original hole has merged in the meantime with another one. If this does not happen 
B' may be called the surface of the same black hole at the later time S. 

Following an idea of Hawking's [24 ], this can be applied to the future evolution 
of multi black hole data as follows. As we already mentioned, the event horizon lies 
on or outside the apparent horizon. Hence the area of the 'surface' (as just defined) 
of a black hole is bounded below by the area of the corresponding apparent horizon, 
which in turn has the lower bound stated in d50t . Suppose that after a long time our 
configuration settles in an approximately stationary state, at least for some interior 
region where no gravitational radiation is emitted anymore. Since our data have 
zero linear and angular momentum, the final state is static and uniquely given by a 
single Schwarzschild hole of some final mass Affln^i and corresponding surface area 
^flnai = 167rAf|^^^,. This is a direct consequence of known black hole uniqueness 
theorems (see |26| or p. 157-186 of [251 for a summary). By the area theorem 
Aflnai is not less than the sum of all initial apparent horizon surface areas. This 
immediately gives 



Y^aj. (51) 



In passing we remark that applied to a single black hole this argument shows that 
it cannot lose its mass below the value rn-i^ := A'°'*''^'/167r, called its irreducible 
mass. Back to the multi-hole case, the total initial mass is given by the straightfor- 
ward generalization of d48l : 



Mnitial = XI = ^ X] ^^^^ + "^^^"^^ X« = X] 1^ I • (52) 
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Using these two equations we can write down a lower bound for the fractional 
energy loss into gravitational radiation: 



AM ^ Mi„itiai - Mfi„,i ^ _ ySi 

For a collision of n initially widely separated {xi 0) holes of equal mass this 
becomes 

AM r-r- 

-^ = 1-^. (54) 

For just two holes this means that at most 29% of their total rest mass can be 
radiated away. But this efficiency can be enhanced if the energy is distributed over 
a larger number of black holes. 

Another way to raise the upper bound for the efficiency is to consider spinning 
black holes. For two holes the maximal value of 50% can be derived by starting 
with two extremal black holes (i.e. of maximal angular momentum: J = m? in 
geometric units) which merge to form a single unspinning black hole l24l . 

One can also envisage a situation where one hole participates in a scattering 
process but does not merge. Rather it gets kicked out of the collision zone and 
settles without spin (for simplicity) in an quasistationary state (for some time) far 
apart. The question is what fraction of energy the area theorem allows it to lose. 
Let this be the kth hole. Then m^"='' > 2afc = m^"'"='V(l + Xk)- Hence 

< < 1 , (55) 

m^— - l + Xk 

showing that an appreciable efficiency can only be obtained if the data are such 
that Xk is not too close to zero. This means that the kth hole was originally not too 
far from the others. This seems an unlikely process. Hence it is difficult to extract 
energy from a single unspinning hole. 

For a single spinning black hole the situation is again different. Spinning it 
down from an extreme state to zero angular momentum sets an upper bound for 
the efficiency from the area law of 29%. This follows easily from the following 
relation between mass, irreducible mass, and angular momentum for a Kerr black 
hole (see e.g. 1.32 J . formula (33.60)): 

= ml + J'^/Aml (56) 



Setting J = m? one solves for m-^^jm = ■y/l/2; hence (m — m;^)/m = 1 — 
1^1/2 0.29. It can moreover be shown fSl that this limit can be (theoretically) 
realised by the Penrose process (compare Frolov's lecture). 
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Figure 4: Multi-Schwarzschild 



Figure 5: Einstein-Rosen manifold 



Needless to say that realistic processes may have far less efficiency than this 
theoretical bounds from the area law alone indicate. Recent numerical studies of 
the head-on (i.e. zero angular momentum) collision of two equal-mass black holes 
give a radiated energy in units of the total energy of only 10~^ |3l. With angular 
momentum the efficiency is of course expected to be much better. Here recent 
numerical investigations give an estimate of 3 x 10^^ for an inspiral of two equal- 
mass non-spinning black holes from the innermost stable circular orbit 12J : still a 
long way from the theoretical upper bound. 

5.5.5 Other topologies 

Other topologies can be found which support initial data with apparent horizons. 
For example, instead of the 'Schwarzschild' manifold with n + 1 ends for n black 
holes one can find one which has just two ends for any number > 2 of holes and 
which has been termed the Einstein-Rosen manifold fl9\. The difference to the 
data already discussed does not primarily lie in the physics they represent. After all 
their different topologies are hidden behind event horizons for the outside observer, 
even though their interaction energies are slightly different 120||22|. However, the 
point we wish to stress here is that such data can analytically and numerically be 
more convenient, despite the fact that the underlying manifold might seem topolog- 
ically more complicated. The reason is that these data have more symmetries, and 
that coordinate systems can be found for which these symmetries take simple ana- 
lytic expressions. For example, in the Einstein-Rosen manifold the upper and lower 
ends are isometrically related by reflections about the minimal 2-spheres in each 
connecting tube, with the fixed-point sets being the apparent horizons. Hence all 
the apparent horizons can analytically be easily located in the multi-hole Einstein- 
Rosen manifold, in contrast to the multi-hole Schwarzschild manifold. 
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Figure 6: bipolar coordinates 



Let us briefly explain this 
for two holes of equal mass. 
Here one starts again from 
coordinatised by spherical bi- 
polar coordinates. These are 
obtained from bi-polar coordi- 
nates {fi, 1]) in the xz plane 
by adding an azimuthal an- 
gle if corresponding to a ro- 
tation about the z axis, just 
like ordinary spherical polar 
coordinates are produced from 
ordinary polar coordinates in 
the xz-plane. The coordinates 

(/U, rf) parametrise the xz-plane according to exp(^ — i-q) = + c)/(^ — c), where 
^ = z + ix and c > is a constant. The lines of constant /i intersect those of 
constant rj orthogonally. Both families consist of circles; those in the first family 
are centered on the z-axis with radii c/ sinh// at z = c coth /i, and those in the 
second family on the x-axis with radii c/| sin i]\ at |x| = c cot rj. 

Following an idea of Misner |31 1, one can borrow the method of images from 
electrostatics (see. e.g. chapter 2.1 in [28 1) to construct solutions ^ to the Laplace 
equation such that the metric hij = ^'^5ij has a number of reflection isometrics 
about two-spheres, one for each hole. In the two hole case one uses the two two- 
spheres fi = ±fiQ for some fiQ > 0, which then become the apparent horizons. 
Using these isometrics we can take two copies of our initial manifold, excise the 
balls 1^1 > fiQ and glue the two remaining parts 'back to back' along the two 
boundaries fi = fiQ and = — ^o- The isometry-property is necessary so that 
the metric continues to be smooth across the seam. This gives an Einstein-Rosen 
manifold with two tubes (or 'bridges', as they are sometimes called) connecting 
two asymptotically flat regions. 

In fact, we could have just taken one copy of the original manifold, excised the 
balls \fi\ > /xq, and mutually glued together the two boundaries fi = ±^o- This 
also gives a smooth metric across the seam and results in a manifold known as 
the Misner wormhole ll30l . Metrically the Misner wormhole is locally isometric 
to the Einstein-Rosen manifold with two tubes (which is its 'double cover'), but 
their topologies obviously differ. This means that for the observer outside the ap- 
parent horizons these two data sets are indistinguishable. This is not quite true for 
the Einstein-Rosen and Schwarzschild data, which are not locally isometric. Even 
without exploring the region inside the horizons (which anyway is rendered impos- 
sible by existing results on topological censorship 1 16 1) they slightly differ in their 
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Figure 7: The Misner Wormhole representing two black holes 



interaction energy and other geometric quantities, like e.g. the tidal deformation of 
the apparent horizons. 

The two parameters c and fiQ now label the two-hole configurations of equal 
mass. (In the Schwarzschild case the two independent parameters were a = ai = 
a2 and ri2 = |xi — 2;2|.) But unlike the Schwarzschild case, we can now give 
not only closed analytic expressions for the total mass M and individual mass m 
in terms of the two parameters, but also for the geodesic distance of the apparent 
horizons I. {£ is used as definition for 'instantaneous distance of the two black 
holes'; for the Misner wormhole, where the two apparent horizons are identified, 
this corresponds to the length of the shortest geodesic winding once around the 
wormhole.) These read 



You might rightly wonder what 'individual mass' should be if there is no internal 
end associated to each black hole where the ADM-formula d35l can be applied. 
The answer is that there are alternative definitions of 'quasi-local-mass' which can 
be applied even without asymptotic ends. The one we used for the expression 
of m above is due to Lindquist f29] and is easy to compute in connection with 
the method of images. But it is lacking a deeper mathematical foundation. An 
alternative which is mathematically better founded is due to Penrose |34|, which 
however is much harder to calculate and only applies to a limited set of situations 
(it agrees with the ADM mass whenever both definitions apply). Amongst them 
are, however, all time-symmetric conformally flat data, and for the data above the 
Penrose mass has fortunately been calculated in l.36j . The expression for m is 




(57) 
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rather complicated and differs from that given above, though the difference is only 
of sixth order in an expansion in (mass/distance) E01 . 

In summary, we see that the problem of setting up initial data for two black 
holes of given individual mass and given separation has no unique answer. Metri- 
cally as well as topologically different data sets can be found which have the same 
right to be called a realization of such a configuration. For holes without associated 
asymptotically flat ends no unambiguous definition for quasi-local mass exists. 

5.6 Non time-symmetric data 

According to a prescription found by Bowen and York Ii37l . we can add linear and 
angular momentum within the setting of maximal data. We can still use confor- 
mally flat data, i.e. set hij = 6ij, on multiply punctured M^. Then the following 
two expressions add linear momentum and spin angular momentum 5* to the 
puncture x = 0: 

k'^ = A (^pin^ + pjn' - {6'^ - n'n^){P ■ n)) , (58) 

K'^ = ^{{Sxfiyn^ + iSxnyn'). (59) 

It is straightforward to check that these expressions satisfy (l38l (note that Z)j = di). 
One can also check that these data will indeed give the proposed momenta and 
angular momenta at infinity (i.e. at the end r oo). For this one may just use 
the 'hatted' quantities in d36t and d37l . since the rescaling d^Tl does not influence 
the leading order parts in the 1/r expansion of K, which alone contribute to these 
integrals. Linearity of all these equations in K allows us to just add Kp and Ks 
and get initial data for one black hole with given momentum P and (spin) angular 
momentum S. Moreover, we can add any finite number of expressions of the kind 
(I58] andl59l with parameters Pi, Si based at the puncture Xi, where z = 1, • • • , n. 
This then leads to a data set whose total linear and angular momentum is given 
by the sum ^ - Pi and Si respectively. But one may not immediately conclude 
that the Pi and Si are linear and angular momenta of the individual black holes. 
Rather, the latter must be calculated for the internal ends of the manifold and for 
this one needs to know <I>. The task then remains to solve ( l39l for the conformal 
factor, with blow-ups being allowed at the given punctures. 

One interesting idea to facilitate solving (l39t is to first split off the singular 
part of <I>, which blows up at the punctures {xi, • • • , x„} like l/\x — Xi\, from the 
regular reminder |6| (compare also |5|). One writes 

1 i v-^ a,- 

^ = - + U, with - := 2^ \ , (60) 

i=n 
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where the > may be freely prescribed. Inserting this into ( l39t gives 

AU + + aUy^ = 0, with (3 = K'^ Kij . (61) 

The point is now the following: for x — > Xj the function a tends to zero as \x — Xi\; 
hence j3, too, tends to zero as \x — Xi\. This means that d^Tl has continuous coef- 
ficients everywhere in since the l/\x — Xi|®-singularity at x-i of the K-squai^ed 
term is cancelled by multiplication with . (Note that this relies on using the /C's 
from d58ll59l . which possess no 1/r" terms with n > 3.) This means that equa- 
tion ( FTl for U can be solved on all of M^, without the need to excise the points 
{xi, • • • ,Xn} and therefore without the need to specify 'inner' boundary condi- 
tions for U; only the 'outer' boundary-condition U{r ^ oo) — > 1 remains. This 
simplification seems particularly useful in numerical implementations (compare 

The total mass of our configuration is M = J2i The individual masses are 
determined just as in l5.5.2l by introducing the inverted radial coordinate fi = af/ ri 
and reading off the coefficient of the 1 / 2fj term in the f j ^ oo expansion. One 
easily gets 

nii = 2ai{U {xi) + Xi) (62) 

with Xi as in d52b . 

The linear and angular momenta at, say, the /cth end can also be calculated by 
using inverted coordinates, given by x = (x — Xfc)a|/r^. Expressed in these co- 
ordinates, the 'hatted' (unphysical) extrinsic curvature tensor is given by JIJ^ Kij 
where := ^ = {a\/f'^)Rl. with R\. = 5\ — 2n*nfc, which is an orthogonal 
matrix. The 'physical' extrinsic curvature is then obtained by multiplication with 
$-2 (compareEJl. Now, ^>(x) = ^(1 + ^ + 0{{l/ff)) so that 

= I [^f + terms oc (^)^} (R'lR^kki), where p>7. (63) 

Inserting the expression ( I58t for Kp results in a falloff so that the individual 
linear momenta are all zero as measured from the internal ends. One may say that 
the asymptotically flat internal ends represent the local rest frames of the black 
holes. Note that these rest frames are inertial since each black hole is freely falling. 
Inserting Kg from ( 15^ gives a falloff and an angular momentum which is just 
—Sk for the fcth end. (Here one uses that R\. is orientation-reversing orthogonal, 
hence changing the sign of eijk, and that R\{S x n)^ = — (5" x n)*.) 
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6 Problems and recent developments 



In this final section we draw attention to some of tiie current problems and devel- 
opments, without claiming completeness. 

I. Given black hole data for n holes of fixed masses and mutual separations 
(whatever definitions one uses here). One would like to minimize these data 
on the amount of outgoing radiation energy. Any excess over the minimal 
amount can be said to be 'already contained' initially. But so far no local (in 
time) criterion is known which quantifies the amount of gravitational radia- 
tion in an initial data set. First hints at the possibility that some (Newman- 
Penrose) conserved quantities could be useful here were discussed in IITbI . 

II. Restricting to spatially conformally flat metrics seems to be too narrow. It 
has been shown that there aie no conformally fiat spatial slices in Kerr space- 
time which are axisymmetric and reduce to slices of constant Schwarzschild 
time in the limit of vanishing angular momentum [18.1 . Accordingly, Bowen- 
York data, even for a single black hole, contain excess gravitational radiation 
due to the relaxation of the individual holes to Kerr form See also fST] 
for an informal discussion of this and related problems. An alternative to the 
Bowen-York data which describe two spinning black holes and which reduce 
to Kerr data if the mass of one hole goes to zero have been discussed in fT2ll . 

III. Even for the most simple 2-hole data (Schwarzschild or Einstein-Rosen) it 
is not known whether the evolving spacetime will have a suitably smooth 
asymptotic structure at future-lightlike infinity (i.e. 'scri-plus'). As a conse- 
quence, we still do not know whether we can give a rigorous mathematical 
meaning to the notion of 'energy loss by gravitational radiation' in this case 
of the simplest head-on collision of two black holes! The difficult analytical 
problems involved are studied in the framework of the so-called 'conformal 
field equations'. See 1 17 1 (in particular section 4) for a summary and refer- 
ences. 

IV. We usually like to ask 'Newtonian' questions, like: given two black holes 
of individual masses mi 2 and mutual separation £, what is their binding 
energy? For such a question to make sense we need good concepts of quasi- 
local mass and distance. But these are ambiguous concepts in GR. Different 
definitions of 'quasi-local mass' and 'distance' amount to differences in the 
calculated binding energies which can be a few 10^^ times total energy at 
closest encounter |20|. This is of the same order of magnitude as the total 
energy lost into gravitational radiation found in |3| for head-on (i.e. zero 
angular momentum) collision of two black holes modelled with Misner data. 
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7 Appendix: Eq. © satisfies tlie energy principle 



By the 'energy principle' we understand the property, that all energy of the self- 
gravitating system serves as source for the gravitational field. In this appendix we 
wish to prove that Q indeed satisfies this principle. For the uniqueness argument 
see l2T]| . 

Given a matter distribution p immersed in its own gravitational potential (j). 
Suppose we redistribute the matter within a bounded region of space by actively 
dragging it along the flow lines of a vector field ^ which vanishes outside some 
bounded region. The rate of change, 6p, of the matter distribution is then deter- 
mined through 6pdV = —L^{pdV) = —V • (^p) dV , where L^-is the Lie deriva- 
tive with respect to ^ and dV is the standard spatial volume element. Note that the 
latter also needs to be differentiated along ^, resulting in L^dV = V- ^. Hence we 

have 6p = — V • (^p). The rate of work done to the system during this process is 

SA= [ dVpfV(t> = -[ dV4>V-{ip)=[ dV(t)6p, (64) 
Jw^ Jm? Jm.^ 

where the integration by parts does not lead to surface terms due to vanishing 
outside a bounded region. Equation (l64t is still completely general, that is, inde- 
pendent of the field equation for (/). The field equation comes in when we assume 
that the process of redistribution is carried out adiabatically, which means that at 
each stage during the process (j) satisfies its field equation with the instantaneous 
matter distribution. Our claim will be proven if under the hypothesis that cp satisfies 
^ we can show that 6 A = c^SMg, where Mq is defined in ^ and represents the 
total gravitating energy according to the field equation. Setting = 4' 

using the more convenient equation we have 



c 



,4 



c4 

lim / dan ■ 



27rG 



r— >oo 



S2(r) 



2ttG 



dV [V'A((5V) - {^^)5i)] 

(65) 



Now, the fall-off condition for r ^ oo imply that Vijj falls off as fast as and 
Sif: as 1/r. Hence the second term in the last line of d65l does not contribute so 
that we may reverse its sign. This leads to 

6 A = 6 lim / da {n ■ = 6 lim / dan ■ \7(f) 

2ttG -r-^^Js^r) 47rG '■^oo752{r) 

= c^6Mg (66) 
which proves the claim. 
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